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Brownian rotors play an important role in biological systems and in future nano-technological 
applications. However the mechanisms determining their dynamics, efficiency and performance 
remain to be characterized. Here the FO portion of the F-ATP synthase is considered as a paradigm 
of a Brownian rotor. In a generic analytical model we analyze the stochastic rotation of FO-like 
motors as a function of the driving free energy difference and of the free energy profile the rotor is 
subjected to. The latter is composed of the rotor interaction with its surroundings, of the free energy 
of chemical transitions, and of the workload. The dynamics and mechanical efficiency of the rotor 
depends on the magnitude of its stochastic motion driven by the free energy energy difference and its 
rectification on the reaction-diffusion path. We analyze which free energy profiles provide maximum 
flow and how their arrangement on the underlying reaction-diffusion path affects rectification and 
- by this - the efficiency. 



Introduction 

Molecular motors play a central role in biological systems where they convert a free energy difference into mechanical 
force. There exist two fundamental mechanisms by which this can be achieved jT]. The free energy difference, coming 
e.g. from a chemical reaction, can be converted directly into mechanical force by some kind of power strike 
The other mechanism is more sophisticated, since here the mechanical movement is provided by thermal fluctuations 
of the motor |SJ El [7]. In the latter case the free energy difference is transformed into entropic forces biasing the 
direction of theses fluctuations . Both mechanism may also be combined and work synergistically jS] . 

An archetype of Brownian rotors is the FO portion of the F-ATP synthase [HI EHl IH] ■ An electro-chemical gradient of 
protons or sodium ions across the inner membrane of mitochondria builds up a free energy difference which is converted 
into a mechanical torque. The latter drives ATP synthesis, the essential energy carrier in living objects. One accepted 
descriptive model for the torque generation of the proton driven F-ATP synthase is briefly (Fig. 1): The FO motor 
consists of a rotary ring, carrying identical hairpin-like protomers in which aspartic or glutamic acid residues transport 
protons along the electro-chemical gradient. The ring is subject to thermal fluctuations and electrostatic interactions 
with its molecular neighborhood. The latter result from the rotor's interface with a positively charged stator, which 
attracts negatively charged (deprotonated) protomers of the rotor. In addition, hydrophilic/hydrophobic interactions 
confine these charged protomers to the hydrophilic stator region [3] [10]. The protomer facing the hydrophobic 
membrane requires release from this constraint through neutralization by a proton coming from the access channels 
[12j . These channels have contact to the respective membrane sides. The frequency by which the constraint is 
removed is proportional to the proton concentrations in the channel, or more precisely to their activity when electrical 
fields are present. A concentration or activity difference of protons between the two channel access sites builds up 
a corresponding probability difference to find one protonated protomer, i.e. occupied proton carrier. This entropic 
force directs proton flow towards equilibration of the proton concentration difference. Rectification of this flow by 
appropriate interactions then produces directed rotation and torque [T?. 

Though there exists valuable and detailed information in literature about the mechanisms of rotation of the FO 
motor |14L 1151 116] and its transduction in ATP synthesis |17j many fundamental questions remain to be solved. How 
does the free energy profile on the reaction-diffusion path, i.e. ring neighborhood interaction, chemical energy and 
workload, mediate the effect of the driving free energy difference on proton flow, and how does it rectify this flow to 
achieve directed rotation? Is there an optimum arrangement and strength of interactions for maximum rotation speed? 
In which way does the binding strength of protons affect the ring rotation? We could recently demonstrate that - 
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O Protonated / deprotonated protomers 
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FIG. 1: Schematic showing the rotation of one protomer unit of the FO portion of the F-ATP synthase, marked with X. Protons 
of either membrane side have access to the ring via channels. The concentration gradient ci > C2 biases counterclockwise 
rotation (see text), resulting in proton flow in the gradient's direction. 



depending on its magnitude - an attractive interaction may support or hamper fluctuation driven motion '18]. What 
is the performance of the ring motor, and when does it reach its maximum? Some of the aspects mentioned above were 
addressed in previous works by simulations or numerical techniques [TTJ [M] in which many details of the biological 
system are considered. In contrast, we will develop here a generic analytical model to obtain fundamental relationships 
that can answer the above mentioned questions. This implies that we consider a model which reduces details of the 
original biological system but still conserves its fundamental properties. Especially for technical applications this will 
be of importance for the development of an optimal working (supra) molecular motor / nano-motor which, triggered 
by chemical breakthroughs [T^, has become an emerging field [IJ. 

The paper is structured as follows: Based on the short qualitative description of the operation mode above, an 
analytic model of the FO motor is presented. For didactical reasons we start with a full rectified rotor, i.e. a perfect 
coupling between proton flow and rotation, which is achieved by the assumption that deprotonated protomers are 
strictly confined to the hydrophiiic stator, while protonated protomers are confined to the hydrophobic membrane 
region [T2^ . This implies a simple circular topology of the reaction-diffusion path, on which the optimum arrangement 
of ring neighborhood interactions providing maximum flow, i.e. rotation speed, is derived. This is done in the 
absence, and presence of external workload. In the latter case also the performance, defined as external work per 
time, is analyzed. In the third section we give up the constraint of perfect rectification, which implies a more complex 
diffusion-reaction path. We analyze the coupling of flow and rotation and the stall workload, determining the efficiency 
of the rotor as a function of the interactions. 



Motor with Check Valve Mechanism 



No Workload 



In our model we consider rotation by one protomer unit as a cyclic process (Fig. 1-2). Two states are distinguished: 
a protonated and a deprotonated one, where the latter refers to the maximum number of negatively charged (deproto- 
nated) protomers within the hydrophiiic region. Within each state the rotor undergoes diffusive motion within some 
free energy landscape (Fig. 2), which is determined by the ring neighborhood interactions, e.g. electrostatic forces 
between deprotonated protomers and the positively charged stator, and hydrophilic/hydrophobic interactions. Tran- 
sitions between the states occur at the locations of the channels, quantified by reaction rates k-^. , k- for protonation 
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FIG. 2: Free energy profile of tlie rotor (deprotonated state red, protonated state green, cliemical transitions black) as a 
function of the cyclic spatial variable x. The location of the channels in which protonation/de-protonation takes place is 
marked by filled circles. Hydrophilic and hydrophobic interactions confine the deprotonated state to the left, the protonated 
to the right section of the diflusion path (see text), and transitions are solely feasible by chemical reaction. This implies that 
the diffusion-reaction path has a circular topology (insert). The free energy course during one cycle on this path is marked 
by arrows {IV ^ I ^ II ■ ■ ■). The choice of proton concentrations ci > C2 makes the rotor run in counterclockwise direction 
and it returns to point IV (hatched circle) at a lower free energy state (blue) — ln(ci/c2). The insert shows the sections of the 
difTusion and chemical transition paths forming a circular topology and the corresponding diffusive flows and chemical fluxes 
Jd, Jd, Ji^u , Jiii^iv ■ 



and de-protonation respectively, 

deprotonated state s ^ protonated state (1) 

where c is the proton concentration in the channel. As suggested in the original works [SJ [TD] and stated in Ref. [T^] , 
we assume in this section that protoniers must be the negatively charged (deprotonated) when facing the positively 
charged hydrophilic stator region, and neutralized (protonated) when facing the hydrophobic membrane . So in our 
model very high energy barriers confine the deprotonated and protonated state to two complementary spatial sections 
which are interconnected by the chemical reaction pathways. This makes the topology of the diffusion-plus-chemical- 
reaction pathway that of an oriented circle (Fig. 2) , and the free energy difference can drive the motor solely in one 
direction of rotation. So the high energy barriers act as mechanical check valves rectifying motion and maintaining a 
perfect coupling of proton flow and rotation. 

Formally the state of the rotor is determined by its protonation state d, p and its local position within this state 
X. The dynamics of the system consists of diffusion and chemical transitions which is described by a set of two 
Smoluchowski equations [20^ to which chemical reaction terms are added. One obtains for the probability density, 

dtPd{x, t) = dxDd{x) [dx - Fd{x)] pd{x, t) - A 

dtpp{x,t) = dMx)[dx~Fp{x)]pp{x,t) + A, (2) 

where Di{x) {i = d,p) is the diffusion coefficient, which we assume to be constant in the following, and the force 
Fi{x) describes the interaction of the rotor with its surrounding. In the one-dimensional case considered here, this 
mechanical force can always be derived from a local potential Fi(x) = —'^[{x). The term A derives from the chemical 
reaction ([T]) and describes protonation and de-protonation via 

A = fc+ ci pd{x, t)S{x - xi^n) - fc_ pp{x, t)S{x - xi^h) + 

k+ C2 Pd{x,t)5{x - xiji^iv) - k_ pp{x,t)S{x - xjjj^jy) , (3) 

where ci, C2 are the proton concentrations in the respective channels. The delta distribution 6{x — xa^b) locates 
chemical reactions on the two cross points of chemical transition and respective diffusion paths which are labeled by 
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the roman subscribes A, B = I, . . . , IV (Fig. 2). Since the rotor must be in some state, conservation of probability 
holds, 

/ pd{x,t)dx+ / pp{x,t)dx = l, (4) 
JxeSd JxeSp 

where Sd Sd are the spatial sections to which the rotor is confined in the deprotonated/ protonated state. 

The stochastic motion of ring rotation is directly related to the diffusive probability flow which in the respective 
protonation states has the form 

j,{x, t) = -D, [d^ - F,{x)] pi{x, t) . (5) 
The chemical flux at the transition between deprotonated and protonated state is 

ji^ii{t) = fc+ ci pd{xi^ii,t) - fc_ pp{xi^ii,t) 

jlll^iv{t) = fc_ Pp{xiii^iv,t) - fc+ C2 Pd{xiii^iv,t) , (6) 

In the steady state the probability densities become stationary pi(x,t) — > pi{x), i — d,p and conservation of flow 
holds. For diffusive flow this implies ji{x,t) = J^. In addition the circular topology of the diffusion-reaction path 
implies that diffusive flow and chemical flux are constant throughout this path, i.e. when we consider flow in direction 
IV ~* I ^ II ^ III ^ IV we obtain 

Jd = Ji^ii = Jp = Jiii^iv = J ■ (7) 

We will now derive the steady state proton flow J, and hence steady state rotation. This requires determi- 
nation of the steady state probability densities at the transition points pi = pd{xi), pu = Pp{xii), pm = 
Pp{xiii), piv = Pd{xiv)- On the diffusion paths the respective gradient of these probability densities pA — Pb 
with {A, B) — {II, III), {IV, I), maintains the corresponding diffusive flow. By generalization of Fick's law for gradi- 
ent driven diffusion we could recently derive this flow compactly as a function of occupation probability in, and first 
passage time through the path [21j. In detail: for a perfect absorbing boundary pB = 0, unidirectional flow in the 
steady state is 

J nA^B 

Ja^B — PA , (,oj 

TA^B 

and vice versa for B ^ A. Here 

TA^B - D-^ n dx e*(") r d^ e-*(«) (9) 

JxA JxA 

is the regular mean first passage time [32], and 



nA-,B = Pa^ / p{x)dx (10) 

JxA 

is the specific occupation number, which by pA tia^b defines the probability to find the system within the diffusion 
path. This number is independent from the boundary value pA, as long as non self interacting diffusing systems are 
considered, which is the case in our model (see Eq. Q. Bidirectional steady state flow for arbitrary values pA, Ps, is 
simply the superposition of unidirectional flows, Ja^b — Ja^b + Jb^a- Since this flow vanishes for equal boundary 
densities |23j, Eq. ^ implies ua^b /ta^b = nB^A/^B^A- The generalized macroscopic Fick's diffusion law then 
reads 



Ja^b = - {pa - Pb) 



(11) 
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with symmetrized first passage time t and specific occupation number n 

1 



^{ta~.b + tb^a) (12) 
^ {nA-,B + ns^A) ■ (13) 



Both parameters depend on the interaction <i> by |18l 

L 



2 
2D 



(14) 
(15) 



where the brackets denote the spatial average () = \/L , with L and D are the length and diffusion coefficient on 
the respective path. 

When we set diffusive flows of Eq. (Ill and chemical fluxes of Eqs. ([6| according to Eq. ([t]) equal J, we obtain a 
set of four linear equations for p/, ■ ■ • , piv- 

Ji^ii = k+ d pi - fc_ pn 
Jp = Jii^iii = np/Tp{pii - Pin) I J ^g-j 

Jiii^iv = fc- Phi - k+ C2 piv ' 
Jd = Jiv^i = nd/Td{piv - Pi) 



In addition conservation of probability (Eq. Q) to find the system within some state must hold, i.e. when expressed 
in terms of specific occupation probabilities (Eq. (lOl) this reads 

Pi ni^iv + pii njj^jjj + pm nm^n + +pjv niy^i = 1 . (17) 

So in summary we have five linear equations from which the four probability densities at the transition points 
P/j''' tPiv E^nd steady state flow J can be obtained. When we assume that the protonation / de-protonation 
reaction rate is much faster than diffusive motion of the motor, i.e. access of protons from the bulk to the rotor is 
not limiting rotation [SF , flow is obtained as 



with 



\K{c^- c) 

Tp + Kci Kc2 Td + K{ci + C2) ^ + K{ci - C2) (^^4^ ' 



At = ^ {ta^b ~ tb^a) , (19) 



as the asymmetric counterpart of the symmetric first passage time r (Eq. 12 1. At quantifies the asymmetry of the 



interaction <i>(x), and vanishes when <i>(a;) is symmetric on the diffusion path. The generalized equilibrium constant 

k+ Up 

^ ^ T. 

K- rid 

/p-*p(2:)\ T 

(e-*<^W) Ld ' ^'^> 

comprises the ring-neighborhood interaction <f> and the standard free energy of protonation 

go = -Hk+/k_). (21) 

Note that the latter is related to the acid dissociation constant of the proton carrier by go « —2.3 pKa- The term 
equilibrium constant becomes evident, when one considers the system under equilibrium conditions, i.e. Ci = C2 = c, 
and hence pi = piv, Pii Piii, and {k+/k^)c = pn/ pi. Then the definition of the specific occupation number 
implies that Kc gives is the ratio of the occupation probabilities N = pn, i.e. Kc = Np{c)/Nd{c). Thus, procedures 
which favor the probability to find the system in the protonated state, e.g. lowering go or an increase/decrease of 



6 




FIG. 3: Proton flow as a function of the constant K (see Eq. ( |20[ )), which comprises ring-neighborhood interaction and proton 
binding strength of the proton shuttle (protomer). The protonated and deprotonated state are assumed to have identical 
diffusion properties {Lp = La, diffusion constant D), and potentials are constant, i.e. Tp = Lp/2D = — Flow is 

normalized to the first passage time ro = {Lp + Ld)^ j {2D), i.e. the time needed to pass the whole diffusion length L = Lp + Ld- 
Two sce narios are shown, ci = 10~®, C2 — 10~* (narrow curve) and Ci = 10~^, C2 = 10~® (wide curve). Note that accordin g to 



Eq. (23 1 maximum flow is located at K ^ = ^C\C2 ~ 10 ^. The dashed curves are the approximations according to Eqs. 1 22 1 
for t: 



le first scenario. 



binding strength in the protonated/deprotonated (<i>p |, 4>dT) state, increase K and vice versa. 



In the hniiting cases of low and high concentrations Eq. ( 18 ) has the form of a Pick's or inverse Pick's diffusion law 
respectively, i.e. 

Jh^gh « ^ ^d"'^"' - Cr') • (22) 

Hence, increasing by increasing the probability to find the rotor in the protonated state increases flow in the 
setting of low proton concentrations, and vice versa for high proton concentrations. This implies an optimum ratio of 
occupation probabilities, reflected by K, for which flow, and hence rotation speed, reaches a maximum (Fig. 3). For 



symmetric interactions, i.e. At = 0, Eq. (18 1 determines this as 

Kmax = \J Tp/Td 1/ ^(ci C2) . (23) 

The above considerations may be formulated in a very elegant way in terms of potentials, which allows generalization 
to the case of external workload. By defining the potentials of the respective driving forces as 

= - HKci) , (24) 

we obtain 



sinh(-AG'/2) 



2 Ta cosh(-AG/2) + ATa sinh(-AG/2) + Tg cosh ((Gi + G2 - ln(Td/rp))/2) 



(25) 



with AG = Gi — G2 = — ln(ci/c2) as the free energy gained during one cycle on the reaction diffusion path (Fig. 2), and 
Xa = [Xp + Xd)/2 as arithmetic, and Xg = ^jXpX^ as geometric mean values of the variable X over the protonated 
and deprotonated states. It should be mentioned that Eq. (25 1 also holds when the standard free energy of protonation 

differs in the channels gj^^ ■ This situation occurs when not only a chemical but an electro-chemical gradient is 
present with Ago as the membrane potential. According to Eqs. ([20 24) this implies two constants if'"'^', K'^^^ and 
corresponding potentials Gi = — ln{K^'^^ q), with the free energy difference AG = — ln(ci/c2) + Ago- 

With the above Equation one can now determine the optimum relation of ring neighborhood interaction quantified 
by <I> and standard free energy of protonation go to achieve maximum fiow, i.e. rotation. In a first step, we vary the 
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FIG. 4: Free energy landscape for a motor working against constant external workload. For simplicity the free energy related 
to internal ring neighborhood interactions $ is not shown. We assume that high energy barriers, acting as check valves, confine 
the motor to the left diffusion path in the deprotonated and to the right in the protonated state. The motor has to afford the 
work Ud ~ ~FLd in the deprotonated and Up = —FLp in the protonated state, i.e. after one complete turn U = Up ~\- Ud, by 
gaining the free energy AG + U = — ln(ci/c2) + U . 



free energy profile without changing the first passage times. This can be achieved by shifting the standard free energy 
of protonation 50, and/or according to Eqs. (9p5l by constant shifts of $p and <^d- Equation (25 1 then predicts flow 
maximum for 



Gi + G2 = ln(Td/Tp) . 



(26) 



The impact of this condition for the free energies Gi becomes best evident when first passage times in the particular 
protonation states are identical, i.e. Gi — — G2. Instead of taking Gi and G2 as opposing driving forces one can 
interpret Gi and — G2 as synergistic driving forces, i.e. maximum flow occurs when driving forces cooperate optimally. 
For non identical first passage times, one has to adjust by ln(Tp/T(i). 



We now vary the first passage times and simultaneously keep the condition (26) fulfilled by appropriate shifts of 



free energy components. The Cauchy-Schwarz inequality then states for the symmetrized first passage time (Eq. (15l) 
that Ti > L1/{2D) {i — p,d), i.e. the first passage time reaches its minimum value wh en t he interaction potentials 
become constant on the respective diffusion paths. Since then Ar vanishes, flow in Eq. (25 1 reaches its maximum 



J m.ax — 



1 



sinh(-AG/2) 



2 TQC0sh(-AG/2) +rg 



(27) 



With External Workload 



The motor is now supposed to works against an external constant force F. As in the previous section, we assume 
that high energy barriers confine deprotonated and protonated state to complementary diffusion paths (see insert in 
Fig. 2). The free energy landscape is that of the previous section with an additional potential term —F x which 
accounts for the external workload. After one cycle the external work is U = Up + Ud = —F{Lp + Ld) (Fig. 4). 



Flow is determined as in the previous section (Eq. 25 1 after adjusting parameters for the external workload (see 
Appendix). Free energies are adjusted to Gi + G2 Gi + G2 + U and the free energy gained after one cycle 
AG^ AG+U = - ln(ci/c2) + U, 



1 sinh 



AG+U\ 



^ — 7 — ^ (28) 



2 cosh {^^) + sinh {-^^) A7, + Tg cosh 



Gi+G2 + U-lniTjr 
2 



We focus on simple constant ring interactions in the respective diffusion paths, which guarantees that the first 
passage times remain constant when interactions are varied (Eq. |9]). The potential G may then be decomposed into a 
component related to the external force and one related to the internal (ring neighborhood) interactions and chemical 
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FIG. 5; Flow under different workload conditions (7 as a function of ring-neighborhood interaction $, and standard free energy 
of protonation go- Proton concentrations are ci = 10^'', C2 = 10^*. The geometric and diffusion parameters are that of 
Fig. 3, i.e. equal for the protonated and deprotonated state. This implies Gext = (Eq. (30l), i.e. maximum flow occurs at 
$p - *d + ffo = l/2(-i7 + ln(ciC2)), i.e. at -16.1 , -17.1, and -18.1, (Eq. ||32|). 



transitions, Gi = Gext + Gi^int, with (see Appendix) 

G^,int = $p-$d + 5o -ln(Q) , (29) 
Ge.t - -ln[(l-e-^-)/(l-e-^-)] . (30) 

This separation makes the interdependence of ring neighborhood interaction (f> and standard free energy go concerning 
their effect on the driving forces Gi^mt evident. Simulations of the FO ring motor working against at constant workload 
by Elston demonstrated for a stronger binding in the deprotonated state, i.e. a decrease of <I>d, that maximum 
rotation rate occurred at a stronger proton binding value, i.e. a decrease of 170. This is directly predicted by the form 
of Gint, since a balanced shift of interaction and standard free energy keeps Gint constant, i.e. in this case the driving 
force at which maximum rotation rate occurred. 

Figure 5 demonstrates that increasing workload U reduces flow and shifts the maximum toward interactions favoring 
the protonated state. This is evident, since maximum flow occurs for 

Gi + G2 + U = ln(rd/Tp) , (31) 

i.e. at 

- *d + 50 = -Gext + ^[-C/ + ln(ciC2) + ln(T<i/Tp)] . (32) 



When flow is considered as a function of workload (Fig. 6), one obtains a more and more parabolic shape of the 
curves, when internal interactions favor the protonated state, i.e. (<i>p — <i>rf -I- 50) decreases. This parabolic behavior 
implies that the motor even under higher workload conditions maintains sufficient performance. For determination of 
this motor performance the following consideration is useful. When a single motor is at some position x, it performs 
the external mechanical work dW — —Fdx after moving to x + dx. For an ensemble of motors, which is in the steady 
state, the number of motors moving dx at position x is Jdt, i.e. the local work performed is Jdt{—F)dx. Performance 
as work per time over the whole length is then 

P =^ I dx {-F) J 
Jo 

= U J . (33) 

Interestingly the motor may run at high performance even when workload increases as long as "tuning" by the 
interaction (<i>p — <i>d + go) is appropriate (insert in Fig. 6). 
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FIG. 6: Flow as a function of workload U for different ring-neighborhood interactions — + go = 

— 16.1 (black), —17.1 (blue),— 18.1 (red) corresponding to maximum flow for U — 0,2,4 (see previous Figure). At work- 
load U — 4.6 — In(ci/c2) = —AG, flow and hence motor action vanishes. The insert shows the corresponding performance 
(Eq. ||33|)with the same color labeling. 




FIG. 7: Free energy profile of the rotor without workload but no check valve mechanism. The labeling is that of Fig. 2. Here, in 
neither protonation state the rotor is confined to sections of the diffusion path by infinite high energy barriers acting as check 
valves. Instead finite energy barriers (magenta in the deprotonated, dark green in the protonated state) only favor sections 
of the diffusion path, the left (red) in the deprotonated, the right (bright green) in the protonated state. This favors the 
rotor to run in counter clockwise direction and this preferred way is marked by arrows. The insert shows the topology of the 
diffusion-reaction pathways and the corresponding diffusive flows and chemical fluxes. 

Motor without Check Valve Mechanism 
No Workload 

We now consider the rotor without the constraint that infinite high energy barriers confine diffusive motion to 
complementary sections of the diffusion path in respective protonation states (Fig. 7). This implies that the rotor 
has the option to move in both directions at the cross points of chemical transitions and diffusion paths. When 
for example the rotor has undergone the chemical transition (protonation) / — > //, it may proceed its way in the 
protonated state either in counterclockwise (right) or clockwise (left) direction. The gain of free energy (— ln(ci/c2)) 
after one completed cycle is the same for both ways. Hence, it is important to stress, that backward rotation, which 
in our example in Fig. 7 is clockwise, is not just a phenomena of fluctuations, but a " thermodynamically allowed" 
way to gain free energy. From the view of the protons this backward rotation can be interpreted as leak flow, which 
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reduces the coupling between overall flow proton flow and directed rotation. 

In the steady state the conservation of flow holds when chemical flux is distributed into diffusive flows on the left 
and right path and vice versa (insert in Fig. 7). In cyclic direction (see also Eq. ([T])) this reads 



Ji^ii — Jp^ + Jp^ — Jiii^iv — + Jd"^ = J 



(34) 



Since the left and right diffusive flows refer to rotation in opposite direction, the effective rotation is 



J ri 



jir) _ j{r) 



(35) 



where the latter relation follows from Eq. ( 34 ) . One can define the fraction / of overall flow J which is directed into 
effective rotation 



(36) 



which can also interpreted as the coupling strength of rotation to proton flow. In the absence of workload the relation 
I/I < 1 holds, since the conservative drift forces Fi{x) = — <i>^(a;) conserve the direction of diffusive flows {Ji^\ j'f'^) 
within a protonation state i = d, p. This direction is then solely determined by the probability gradient i.e. in our 



setup with ci > ca it points from // III and IV I. Hence, \Jy - J^'\ < max(|4 |, \J^'\) < \ J\. The 
relation |/| < 1 does not hold generally for non-conservative forces, e.g. when a constant external workload is present. 
Depending on its strength this force may rotate the motor within the protonated state, and by this decouple rotation 
Jrot from proton flow J, a situation we will discuss below. 



Since in the steady state the generalized Pick's diffusion law (111 holds for the particular diffusive flow components 
on the left or right paths, 



j(l,r) _ {l,r) i(l,r), n \ - J 



(37) 



with {A,B) = {IVJ), (11,111) and i = p, d, one obtains with Eqs. (14 15 1 



/ = 



L 



(38) 



where L^^\ L^'') 
particular paths. 



I.e. 



are the lengths of the left and right diffusion paths, and () denotes the spatial average on the 
There is no rotation / = 0, when free energy landscapes of the protonation states are congruent, 



)) 



(e*^")/(4' 



Conversely free energy profiles favoring complementary paths 



in respective states synergistically increase rotatory efficiency |/| 1, e.g. as in Fig. (7) a barrier on the right 

(r) (I) 

of the deprotonated state t, and left of the protonated state t- Infinite high barriers imply maximum 
efficiency f — I since they act as check valves, a situation discussed in the previous section. In the real motor (Fig. 1) 
the synergism of free energy profiles is realized by the high hydrophilic-hydrophobic energy barrier which confines 
deprotonated protomers to the hydrophilic stator region. Since there is at least one deprotonated protomer in either 
protonation state, backward rotation is impeded, i.e. f — 1- However, the situation would be different for motors in 
which the maximum number of deprotonated protomers facing the hydrophilic stator is one. In this case there would 
be no mechanism preventing the rotor from backward rotation in the protonated state. 

The steady state flow J can be derived as in the previous section, based on conservation of flow (Eq. ([t]) and overall 
probability (Eq. (17). To apply the Eqs. (16) for this derivation, one has to formulate diffusive flows in the form of 
a generalized Fick's diffusion law Eq. (111. This is not restricted to diffusion on a linear path. Instead it holds for 
steady state diffusion through arbitrary domains, as long as the forces acting inside are conservative, i.e. when they 
derive from a potential F{x) — —^'(x) |21j . Here the domain is no longer a simple path as in the previous sections, 
instead it consists of a left and right path in either protonation state as the insert in Fig. 7 shows: red/magenta 
in the deprotonated, green/dark green in the protonated state. The boundaries of the domains are the cross points 
{A, B) = (/V, /), (//, ///). The application of Eq. ( 11 ) requires the knowledge of the specific occupation number and 
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FIG. 8: Left; free energy topology for a motor working against constant external workload. A high energy barrier confines the 
motor to the left diffusion path in the deprotonated state. In this state, the motor first has to afford the work Ud {IV I) 
before free energy is gained by protonation (7 — > //) . However in contrast to Fig. 5 the motors has in the protonated state the 
choice either to complete the cycle in counterclockwise direction and to afford the work or it may rotate backwards and 

gain back the energy Ud = Up \ So the bifurcation of the reaction-diffusion path at point II implies the option of two paths, 
by which the motor may reach its initial position {IV') and gain free energy. For simplicity the internal interactions <I>p are 
not shown on the left, but on the insert right which has to be superimposed on the free energy profile on the left. 



first passage time. The specific occupation number is simply the sum of the specific occupation number of the left 
and right diffusion path. In both protonation state follows from Eq. (14 1 



,(0 



L('V2 (exp(-$('))) +L('^V2 (exp(-$(''))) 
+ 2,M)/2 (exp(-$)) 



(39) 



To obtain the first passage times one h as t o keep in mind that Eq. \\\\ holds for the particular fiow components on 

either left or right diffusion path (Eq. (37 1). Since Jd^p = Ja^b — ^a^b + '^aI^b^ ^'^^ ^ ~ + "^'^^ inverse 
first passage time derives as the occupation probability weighted sum of the inverse first passage times on the left anf 
right diffusion path 



1 



1 



(40) 



Insertion of specific occupation number and first passage time from the above Equations in Eqs. (11 16 1, then provides 



as in the previous section fiow J from Eqs. (18 25 1 



With Workload 



Now the rotor is now supposed to work against a constant force. For didactical reasons and to avoid too much 
complexity without gain of physical information we will not consider the topology of the reaction -diffusion path in 
its most general form as in Fig. 7. Instead we will confine the rotor by high energy barriers to the left path in the 
deprotonated state. In the protonated state it has the choice to take both ways (Fig. 8). In the presence of an external 
force this option has strong implications on the topology of the free energy landscape. When check valves were present 
this landscape would be just a simple, non-single valued function of a 1-D cyclic reaction diffusion path, i.e. after 
one turn the motor gains the free energy AG -I- U (Fig. 4). However, now the motor has the option to diffuse either 
along the left or right path in the protonated state, i.e. it may gain the free energy AG -I- J7 on the right vs. AG on 
the left path (Fig. 8). In other words this bifurcation of the reaction diffusion path leaves the option for the motor 
to proceed its way in clock- or in counterclockwise direction, however external work is only performed in the latter 
case. In addition gain of free energy is even higher when taking the "ineflicient" path AG < AG + IJ . The effective 
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FIG. 9: Flow, which is effectively transformed into counterclockwise rotation, i.e. against the external force, as a function of 
the standard free energy of protonation go = — ln(A;+/fc_). Concentration are ci = 10~^ and C2 — 10~*, and potentials on 
the paths forming the way in counterclockwise direction are $^ = = 0. The workload is (7 = 2. A high energy barrier 
"I>; = 15 prevents the motor from backward rotation (black line), whereas no barrier = (blue line) implies clockwise rotation 
Jrot < when the energy level of the protonated state is lowered {go decreases). 



^^.0,^0 1.2 




FIG. 10: Flow, which is effectively transformed into counterclockwise rotation Jrot as a function of workload. Concentrations 
and potentials are as in the previous Figure. The standardized free energy of the protonation was chosen to be go — —16.1. 
With an energy barrier favoring counterclockwise direction ($p' = 15), the whole concentration gradient related free energy 



|AG| — ln(ci/c2) — 4.6 may be converted into work. When there is no energy barrier ("l>p' 
as the zero point of this function (Eq. (42|) demonstrates. 



0) the efficiency is much lower. 



rotational component as determined from Eq. (351 is 



J, 



rot — 'Jd '^p 

- 7„ 7(0 - /M- 

— J Up — Up , 



(41) 



j: 



where we exploited that the motor is confined to the left difi^usion path in the deprotonated state, i.e. J — 
The efficiency factor / 



{I) 



(1 — Jp^ / J) has not the simple form as in Eq. (38 1, since the external workload produces 



a non-conservative force field in the protonated state. For the determination of fiow J and its effective rotational 
component Jrot see Appendix. 

We first consider the case that there is no energy barrier in the protonated state which could bias the direction of 
rotation, i.e. ^p"* — (Fig. 9). Decreasing the standard free energy makes effective rotation first pass through a 
maximum and then vanish. At this point proton and leak fiow are equivalent J = JjP and fiow on the right path in 
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ln( c, / 




11: 



Workload Uo for which rotation vanishes as a function of the barrier $p 

6 in-8 



m 



according to Eq. (|42[). Standard free energy 



FIG. 

is go = —16.1 and concentrations are ci = 10~"^, C2 = 10"**. High barriers, which favor counterclockwise rotation against the 
external force {$p^ — + +oo), make the workload approach the free energy difference Uo —^ |AG| = ln(ci/c2) = 4.6. 



(r) 

the protonated state is absent Jp — (Eq. (41 1. Further lowering of go decreases the transition probability from the 
protonated to the deprotonated state and the rotor prefers to remain in the protonated state. The rotor runs along 
the direction of the external force, J^ot < 0; and by moving along the circular free energy landscape of the protonated 
state it gains the free energy —U per round. We assume now constant interaction potentials on the left ^p'' and right 
4>p section of the diffusion path in the protonated state, and in the deprotonated state ^d- At the threshold value 
go for which Jrot vanishes, workload U and ring neighborhood interactions satisfy (see Appendix) 



l + e*p 



= ' (,) . (42) 

This Equation also makes the interdependence between standard free energy and the barrier evident. The 
higher the barrier, the lower is the threshold for go to make rotation vanish, i.e. a high barrier counteracts the external 

(r) 

force. Interestingly the interaction on the right site in the protonated state <I>p does not appear in the equation, 
i.e. it has no effect on the direction of rotation. One could argue that inversely to the barrier effect of a high 
should facilitate negative (clockwise) rotation in direction of the external force. For analysis let us first consider 



a situation in which the rotor runs against (counterclockwise) the external force, i.e. Jrot > 0. Equation (41 1 then 



implies jjf ^ > 0. When we introduce a barrier on the right site by elevating we reduce Jp^^ but it does not 
vanish, or become negative, since this barrier also works against the clockwis e ro tation in direction of the external 

(r) (r) / J l\ 

force. In the limiting case <i>p — > +oo, Jp , and hence Jrot approach zero (Eq. (41 1), but they never become negative. 



(r) 

i.e. there is no change of direction of rotation. So in summary <i>p has no effect on the direction of rotation. The 
above Equation also defines the maximum workload Uq for which rotation against the external force is maintained. 
A workload above this threshold implies backward rotation. When the barrier height ^p^ is raised, this threshold 
increases, and for ^p^ +oo it asymptocially reaches the free energy difference Uq —AG = ln(ci/c2) (Fig. 10). 



Discussion and Conclusion 



Whereas past interest of theorists mainly focused on the working principles of linear motors |28l IHU] , in this 
paper we presented a comprehensive analytical model of FO-like Brownian rotors, which are archetypes of rotatory 
biological molecular motors. These rotors transform a free energy difference maintained by an (electro)-chemical 
gradient into mechanical work by means of random Brownian motion. In contrast to the classical flashing ratchet 
concept |26j . in which on a molecular level chemical energy of the motor is directly transduced via a power strike 
mechanism into mechanical energy jH [3T] , these types of motors convert the free energy difference into entropic forces 
[5J [271 131]. This is realized by (chemical) transitions between different free energy profiles, the latter biasing motion 
into certain directions. Directed average motion results from different transition rates at different locations. In the 
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case of the FO rotor this is realized by different cation (proton) concentrations in respective access channels. Note 
that these transitions - though they may alter the free energy of the motor state - do not alter the mechanical energy 
at the point of transition [32] • In contrast, the latter would be the basic mechanism of a classical flashing ratchet. 

To make a Brownian rotor run optimally requires both: The free energy difference must be converted into maximum 
average motion, and - in addition - this average motion must be rectified into one direction, as in linear motors [33]. 
We demonstrated that in our system free energy may be gained by forward as well as backward rotation. This is 
important since it demonstrates that backward movement is not only due to fluctuations. In our case it exists also as 
average motion on certain segments of the diffusion path - ie. the leak part - and is a consequence of the topology of 
the diffusion-reaction path and the overlying free energy profile. 

Theses considerations stress the importance of the interaction of the rotor with its surrounding. For FO-like rotors 
we derived analytical expressions for the dependence of flow on interactions, determining maximum flow for a given free 
energy difference. It could be demonstrated that even under high workload conditions optimum interactions implied 
a high performance of the motor. Furthermore the arrangement of interactions with respect to its rectifying property 
was analyzed. The latter determines the maximum workload the rotor can accomodate. However, only under perfect 
" check valve " conditions the maximum workload approaches the driving free energy difference. In its biological 
realization this situation is almost accomplished, since in either protonation state at least one deprotonated protomer 
is confined to the hydrophilic stator region (see Fig. 1). This problem of rectifying is typical for molecular rotors 
that work in the strong friction limit, i.e. when motion is over-damped and inertia effects are absent. Macroscopic 
rotors may overcome this problem by exploiting in addition the rectifying effect of inertia, e.g. by some fly wheel 
mechanism. 
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APPENDIX 



In our model the motor gains free energy by running through an assembly of diffusion and chemical transition paths 
(Figs. 2, 4, 7, 8). In the absence of external forces diffusion connected cross points of diffusion and chemical transition 
{IV, I); {II, III) exhibit the same free energy level which allows to write diffusive flow in the form of a Pick's diffusion 
law Eq. (jni This was the base to derive steady state flow throughout the system, i.e. diffusive flow and chemical 
flux (Eqs. (18] 25|. When the motor has to perform a constant workload, i.e. a constant external force is present, 
these cross points exhibit different free energy levels, which in the general case depend on the path the motor takes 
(Fig. 8). To integrate this in a more general derivation of flow one has to start with the basic Equation (|8|, which 
relates unidirectional diffusive flow in the steady state through some domain to specific occupation number and first 
passage time. Unidirectional means that the domain has two boundaries, one acting as a particle source with constant 
concentration/probability density, the other as a pure absorber. In our model the domain consists of the diffusion 
paths connecting two cross points, which act as boundaries. Since no self interaction is assumed, bidirectional flow 
between the boundaries is just the superposition of respective unidirectional flows, i.e. 



Ja^b = PA Pb 



ta^b 



tb- 



(A-1) 



with {A,B) ~ {IV, I), {1 1, 1 1 1). The above Equation is a generalization of Eq. (11) and is valid in general for 
all type of drift forces, i.e. conservative and non-conservative, which act within the domain |21j . In the steady 
state flow, i.e. chemical flux Jj^u, J m^iy (Eq. ([6|) and diffusive flow Jjv^i — Jd, Jii^iii = Jp is constant 
throughout J (Eq. ([t])). Hence with Eq. ( A-1 ) for the diffusive flows and Eqs. (|6| for the chemical fluxes and Eq. ( 17) 



for the conservation of probability we have a system of five linear equations determining the the four probability 
densities at the cross points pi ■ ■ ■ piv and flow J . When we assume that the chemical reactions, i.e. protonation and 
de-protonation, are much faster than rotor diffusion, flow is obtained as 



J = 



Cl 



nii^iii niv^i 



C2 



Tl^IV TlII^II 



Cl I niv^i H nii^iii + C2 nj^iv H nm^u 



TlI^III 



V Tin- 



k- f nj^iv , niv^i \ , k+ ( nu^m nm^u 

1— njv^i H nj^iv + CiC2 -;— njji^jj H njj^jjj 



(A-2) 



The above Equation is valid for all diffusion path topologies connecting the cross points {IV, I) or (//, ///). 



Rotor with Check Valve Mechanism 



We will first consider the more simple case that high barriers act as mechanical check valves and confine the rotor 
in the protonated state to the right, in the deprotonated state to the left path (Figs. 2, 3,). This arrangement avoids 
that the free energy becomes a multi-valued function of the diffusion path in either protonation state, a case which 
we will discuss later on. Hence, the free energy level (p of diffusion connected cross points differs by the workload one 
has to afford to get from one point to the other. With constant force F on the respective diffusion paths one gets 
= Ud = —FLd in the deprotonated, and 0/// — 0// — Up — —FLp in the protonated state (Fig. 4). Diffusive 



flow vanishes when the probability densities exhibit the same chemical activity, i.e. piy 
Together with Eq. (A-1) this implies that 



e'^'^pi and pn 



e^'Piii- 



nA^B 
ta^b 



tb^a 



(A-3) 
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where {A,B) = {IV, I) or {II, III). This suggests to normahze the specific occupation number by the free energy 
level, fiA^B = e^'^'^nA^B- With symmetrized quantities 



T = - {ta^B + Tb^a) 



(A-4) 
(A-5) 



one obtains from Eq. (A- 3) 



riA^B _ nB^A _ n 
TA^B tb^a t 



(A-6) 



This enables us to formulate a Pick's diffusion law, which is similar to that of Eq. (Ill, except that probability 
densities are replaced by corresponding activities 



Ja^b - - {e'^^pA-e'^^pB) 



(A-7) 



Insertion of the normalized specific occupation numbers in Eq. (A-2|, and considering Eqs. ( A-6j|A-7 1 leads to 



(A- 



Here U — Up + Ud denotes the workload performed after one complete cycle, and An — l/2{hA^B — tt-b^a)- One 
can express the above relation in terms of potentials. The driving forces are 



G,; = -In 



= 90- In 



- ln(c,) 



(A-9) 



where go is the standard free energy of the protonation process. Then Eq. (A-8l reads 



sinh - {AG+U)/2 



Ta cosh((AG + U) /2) + sinh(-(AG + U) /2) Ar^ + Tg cosh 



Gi+G2 + i7-ln(rrf/rp) 



(A-10) 



Here Xa = l/2{Xd + Xp) denotes the arithmetic, Xg = ^XdXp the geometric mean of some parameter X in 
the deprotonated and protonated state, Ar =1/2 {ta-^b — tb^a) is the antisymmetric counterpart to the sym- 
metrized first passage time in Eq. (A-5 1. In the derivation from Eq. (A-10 1 from Eq. (A-8l we exploited the relation 
(Arii/ni) = Ar^ for each protonation state i = d p fM\ . 



Specific Occupation Number and First Passage Time in a Linear Potential 



A constant external force F superimposes a linear potential U {x) = —Fx on an existing internal potential, de- 
termined by the ring neighborhood interaction. At first we will now focus solely on the first, and neglect internal 
interactions. We consider unidirectional steady state diffusion on a path in the protonated or deprotonated state 
characterized by the index i = p, d, against the external force from crosspoint A at position x = to crosspoint B at 
position X — Li. The occupation number ha^b gives the probability to find the system within that path normalized 
by the boundary probability density pA- From Eq. ([5| one can readily determine the diffusive flow in the steady state 
by integration when one respects the equivalence {u^ — F{x)) = exp{—U{x))dxexp{U{x)). Within a linear potential 
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one gets 



D 



D 



dx exp(J7(a;) 



e^(°) PA 



PA 



where Ui = U{Li) 



-FLi. The probabiHty density in the steady state is obtained from Eq. ([s]) as 



(A-11) 



p{x) 



_ Ja^b Li 



(e^(-) ~ l) 



(A-12) 



Insertion of the flow of Eq. (A-11 1 and integration of the probabihty density then provides the specific occupation 
number as 



nA^B 



Pa 



dx p{x) 
1 



(A-13) 



Note that for the normalized occupation number ua^b = n-A^B holds, since f/(0) = 0. Vice versa one obtains ub—^a 
by just inverting the sign of Ui in the above Equation, however for the specific occupation number one has to consider 



that fiB^A = e 
state 



ub^a- Hence, one obtains for the symmetrized specific occupation number in either protonation 



2 {flA^B + flB^A) 

I'm 



(A-14) 



Similarly the antisymmetric part is determined as 

An,; = 



2 



nB^A) 



Li sinh(L/",) - Ui 
a e^. - 1 



(A-15) 



Up to now we only considered the influence of the external interaction on the occupation number. Potentials reflecting 
internal interactions have to be superimposed on U{x). In the manuscript they were assumed to be constant ^i on 
the path < a; < L^. Then, the occupation numbers scales with the factor 



ni{<i>i) 



n,(0) 



(A-16) 



as substitution of U{x) U{x) + ^i in Eqs. ( A-11|[A^T3 I demonstrates. The same is true for An. Superimposing 
constant potentials on the diffusion pathways allows factorizing of the generalized equilibrium constant into a com- 
ponent characterizing the internal and external interaction. Therefore the corresponding potential (Eq. A-9 1 may be 
written as 



= 50 - In 



nd 



- Hcj) 
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-In 



1 - exp(-t/p) 
1 - exp{-Ud) 



where we exploited that Up/Lp = Ud/Ld — —F. 

The regular first passage time to pass an interval of length Li from one end A to the other B is [22] 



(A-17) 



dri e 



di e 



where the diffusion coefficient D was taken as constant. For a linear potential one gets 

L? e^^ -[/,;- 1 



TA^B 



D 



Uf 



(A-18) 



(A-19) 



and vice versa for tb^a by changing the sign of Ui. Hence, the symmetrized = 1/2{ta^b + tb^a) and anti- 
symmetrized At^ = 1/2{ta^b ^ '''b^a) fii'st passage times in either protonation state are 



At,: = 



L[ smh\U,/2) 
2D (f/,/2)2 • 
Lf sinh(£/i) - £/, 



(A-20) 



Note that the first passage time in Eq. (A-18 1 does not change when constant potentials <i>i are superimposed, since 
^u(n)+'^i g-c^(4)-*i = e^^"^") e~'^^^\ So the symmetrized and anti-symmetrized first passage times are not affected 
either. 



Rotor and Effective Rotation - General Case 



In the case we discuss in the main paper, the free energy in the protonated state becomes a multi-valued function 
of the diffusion path. When we consider in Fig. 8 diffusion in the protonated state from // to ///, the motor reaches 
the potential 



'^III 



Jjir) 
P 



(A-21) 



when it takes the right path, but 



'Ud 



(A-22) 



when it takes the left path. Hence, one cannot simply as sign cross points free energy levels so that diffusive flow i s 



driven by an activity gradient {e'''" pn — e'''"' pm) (Eq. A-7l a prerequisite to obtain flow from Eqs. (A-8 A-lOl 
However, one can accomplish this by introducing effective potentials 
protonated state is the sum of flow on the left and right path 



(e//) ^ Q^g ]-^j^g keep in mind that flow in the 



J, 



ii^iii 



4' 



j(r) 



(A-23) 



Flow on each particular path is obtained from Fick's Equation (A-7l 



Op 



{exp{(l)fj)pn - exp{(j)fjj)piii) 



Up' 



(A-24) 
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and 



Up 

Tp 
~(r) 

rVp 



{pii - exp(J7^''))p7//) , 



(A-25) 



where we defined the cross point (/J) as reference point with 



0, i.e. only the free energy at point (///) 



becomes path dependent. The parameters fip "^ and Tp ^ are obtained from Eqs. (A-5l. When we define 



= hi 



= hi 



4" 


exp((; 










~ (i) 


'p 


4" 






e p 











(A-26) 



(note U^^ = C/d) we obtain with Eqs. ( A-23|A-24 | 



III 



exp(4y^^) PI I - exp{<l)f//') PI 1 1 



(A-27) 



i.e. difi^usive flow is driven by a gradient of effective activities. It is worth to note that detailed balance is not fulfilled 
in the protonated state, since the external force superimposes a non-conservative force field. A vanishing diffusive 
flow in Eq. (A-271 for equal effective activities just implies that Jp^ = —J^^ i.e. a circular diffusive flow. This will 
be of importance for determination of the effective rotation. 

To bring Eq. (A-271 into the form of Eq. (A-7l, one has to define appropriate expression for the specific occupation 

as 

(A-28) 



number n and first passage time r. According to Eq.(A-l I diffusive flow in the protonated state may written as 

niii-^ii 
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Similarly as in Eqs. ( A-3|A-6" |, a vanishing flow for equal activities (Eq. (A-271 implies from Eq. (A-28 1 

TlI^III 



TlII^II 

Hence, with symmetrized normalized speciflc occupation number and flrst passage time 

1 
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II jnii^iii 



exp( 

(tii^III + Tiii^Il) , 
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t>iii jniii^ii 



flow can be written finally in the form of Eq. ( A-7 1 
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(A-31) 



What is left is the determination of the specific occupation numbers nii^m niii^ii and first passage times 
Tii^iii Tiii^ii. Specific occupation numbers in the protonated state are the sum of respective numbers on each 
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path, nji^jii = n^jj^jjj + 'iT'i'i^mi Emd similarly for njjj^jj. For determination of the first passage times one 
exploits that flow is sum of flow on the left and right path, i.e. for // — > /// 



nii^iii 

TlI^III 



'-ii^iii 
'ii^iii 
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'^II^III 
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and, hence, I/th^ih is just the occupation number weighted average of first passage rates on the single paths 

^/tu^iii- The same hold for parameters in direction /// —> II. The occupation number and first passage times 
characterizing the left or right path are determined according to the previous section. 

Equation (A-311 expresses diffusive flow on the assembly of left and right diffusion paths in the form of Pick's 
diffusion law. Hence, one can directly apply Eqs. (A-8 |A-10| to determine flow in the steady state. Note that the 
workload U in Eq. (A-lOl is replaced by an effective workload 



(A-33) 



The efficiency factor / determining the fraction of flow J transformed into effective rotation, takes, respecting the 
constraint that the motor is conflned to the left path in the deprotonated state ( J^''' = J), the form 



/ = (J("l- J«)/J 
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The fact that different activities act as driving forces on the left and on the left plus right path in the protonated state 
explains that / has not the simple form as in the case without workload (Eq. pS])). Instead one has to determine the 
steady state probability densities at the cross points of chemical transition and diffusion px, X = I, - ■ ■ IV, which 
determines / as 
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where the driving forces Gi are defined in Eq. (A-9), i.e. AG — — ln(ci/c2). Note that in the absence of external 
forces [U^'^^^'^ = 0) the last summand in the bracket vanishes and / takes the simple form of Eq. (38 1. 

When we assume constant interaction potentials in the deprotonated '^d and on the left and right path in the 
protonated state , and insert parameters for work against a constant force (see previous section), we obtain 

/ as a function of the workload the motor performs after one complete cycle U = Ud + Up^ (see Fig. 8) , 



e*r+f^"effo(c2 - ci)s4i|^ + {e^ - eU»){e'^'^ + effo-i"(ci)+*<" ) _ e*^(eC^ - 1) 

e"p -1 
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2sinh(^^)e* , 
^ 2 ' (A-36) 



e4^'e.o(c2 - cO^^g^ + 2 sinh (5^-^) e* - 2 e*^ sinh([//2) 
where 

+ ln(e*''+e*^"+9°-^"^i)) , (A-37) 
and Uq is the external workload when rotation vanishes (/ = 0), 

f P^d I p4>^'^+go-lnc2\ 

C/o = In • (A-38) 

Note that stabilizing the protonated state by qq — > -co or $d +(X) implies /7o ^ 0, $ ^ $d and, hence, 
/ —00, i.e. the rotor spins in the protonated state solely driven by the external force. Stabilizing the right 
path of the protonated state <I>p^ —^ +oo implies that the whole chemical gradient is used for directed rotation 
Uo ln(ci/c2) = -AG and / -> 1. 



